Casimir interaction between a perfect conductor and graphene described by the Dirac 

model 
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We adopt the Dirac model for graphene and calculate the Casimir interaction energy between a 
plane suspended graphene sample and a parallel plane perfect conductor. This is done in two ways. 
First, we use the Quantum Field Theory (QFT) approach and evaluate the leading order diagram 
in a theory with 2+1 dimensional fermions interacting with 3 + 1 dimensional photons. Next, 
we consider an effective theory for the electromagnetic field with matching conditions induced by 
quantum quasi-particles in graphene. The first approach turns out to be the leading order in the 
coupling constant of the second one. The Casimir interaction for this system appears to be rather 
weak. It exhibits a strong dependence on the mass of the quasi-particles in graphene. 

PACS numbers: 12.20.Ds, 73.22.-f 
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I. INTRODUCTION 

Graphene is a (quasi) two dimensional hexagonal lat- 
tice of carbon atoms. At present, it belongs to the most 
interesting materials in solid state physics in view of its 
exceptional properties and its potential applications in 
nano technology (see reviews, Refsji^). At small separa- 
tions, nearly down to contact, the interaction between a 
graphene sample and any solid body (dielectric, conduc- 
tor or another graphene) is due to the van der Waals and, 
at larger separations, Casimir forces (see, e.g. Refs4^). 
The latter are the subject of the present paper. We con- 
sider the same geometry as in the original Casimir effect 
(2 parallel planes) with one plane being graphene and 
the other one - ideal conductor. This setup was consid- 
ered in Refs£>&£ using a hydrodynamical model for the 
electrons in graphene following Refs. 8 ' 9 . Later it became 
clear that this model does not describe the electronic 
properties specific to this novel material. Here we use a 
realistic and well-tested model where the quasi-particles 
in graphene are considered to be fermions subject to the 
Dirac equation which models their linear dispersion law. 
Details of this model can be found in Refsji^. 

The result presented below is the first calculation of 
the Casimir interaction of graphene made within a rea- 
sonable theoretical model. Although the hydrodynamical 
(plasma) model is not applicable to graphene, it works 
well for some other materials, and it will be used as a 
theoretical reference point to compare our results. 

We like to mention some related works. The Casimir- 
like interaction between adatoms due to fermionic modes 
inside graphene was studied in Ref^i. The graphene- 
metal interaction at separations 2-4A (much smaller than 
the Casimir distances) was investigated in Refs^. 

Let us formulate the model. The electronic properties 
of graphene are well described by the above mentioned 



Dirac models, ft incorporates the most essential and 
well-established properties of the quasi-particles' spec- 
trum: the linearity, a very small mass gap (if any), and 
a characteristic propagation velocity which is 300 times 
smaller than the speed of light. The model deals, there- 
fore, with light fermions in 2 + 1 dimensions (confined to 
the surface of graphene) with the action 



d ' xip(j (idi - eAi) - m)V>, 



(1) 



where I 



7°, 7 1 ' 2 



0,1,2. The matrices j l are rescaled, 7 = 
vf! 1 ' 2 : 7o = -(7 1 ) 2 = ~(7 2 ) 2 = I- V F is the 
Fermi velocity. In our units, % = c = 1, vf — (300) -1 . 
The gamma matrices are taken in the form of a direct 
sum of two incquivalent representations (differing by an 
overall sign). There is an additional ('valley') degeneracy 
in graphene, so that we have 4 two-component spinors or 
2 four-component spinors in graphene. The value of the 
mass gap parameter m and mechanisms of its generation 
are under discussio n 13 ' 14 i 15 i 16 . The upper limit on m is 
about O.leV, but probably is much smaller. Due to this 
smallness of the mass, the quasi-particles exhibit a rela- 
tivistic behavior at rather small energies, which makes 
QFT a more adequate language to describe graphene 
than Quantum Mechanics. As we shall see below, the 
Casimir force is very sensitive to the value of m, and thus 
the Casimir experiments may be used to estimate it. 

The propagation of photons in the ambient 3 + 1 di- 
mensional space is described by the Maxwell action 



S 



1 



M 



d A xF^F^, n, u = 0,1,2, 3. (2) 



The coupling constant is normalized according to 
e 2 /(47r) = a ~ 1/137. The Dirac model with quantized 
fermionic quasi-particle excitations and classical electro- 



2 



magnetic modes describes rather well the optical prop- 
erties of graphene. By construction, this model should 
work below the energy scale of about leV, but even above 
this limit the absorption of light by suspended layers of 
graphene is reproduced with a high precision^. 

In the following we calculate the Casimir force be- 
tween a flat suspended monolayer graphene sample and 
a parallel flat perfect conductor. We shall suppose that 
the graphene sample occupies the plane x 3 — a > 0, and 
the conductor corresponds to x 3 = 0. 



II. QFT APPROACH 

One possible way to calculate the Casimir energy in the 
system in question is to evaluate the effective action V in 
a quantum field theory described by the classical action 
Sb + Sm ■ Since the background is static (the positions of 
the surfaces do not depend on time), the energy density 
per unit area of the surfaces is £ = —F/(TS), where, 
because of the translation invariance, one has to divide 
the effective action T by the (temporarily introduced) 
time interval T and the area of the surface S. At the 
leading order in the fine structure constant a we have 



(3) 



where the solid line denotes the fermion propagator in 

2 + 1 dimensions (i.e., inside the graphene sample), and 
the wavy line is the photon propagator in the ambient 

3 + 1 dimensional space subject to the perfect conductor 
boundary conditions 



Ao\ x 3=o — Ai\ x 3 =0 — A 2 \ x 3 =0 — d 3 A 3 



= 0. 



(4) 



We use the Feynman gauge such that no contribution 
from ghosts appears. For details on diagrammatic nota- 
tions see Rcf.— . 

The calculation of ([3]) is rather similar to that of the 
radiative corrections to the Casimir energy in Ref.— . The 
important difference to those works is that they consid- 
ered the fermions propagating in the ambient space as a 
correction to the Casimir force between two perfect con- 
ductors, which appeared to be very small. In that case 
the diagram ([3J represents only the 0(a) correction to 
the Casimir energy. In our case, the diagram ([3]), still 
being 0(a), represents the leading order effect. 

The fermion loop in 2 + 1 dimensions has already been 
calculated in a number of paper s 13 ! 15 ' 16 . This result is 
most conveniently expressed in terms of the effective ac- 
tion for fermions in the presence of an external electro- 
magnetic field, S e s(A) = —i\ndet(('j l (idi — eAi) — m). 
To the quadratic order in A the effective action reads 

S cS (A) = A w^^w A 



d 3 p 



A J (p)W l (p)M P ), 



where 
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a$(p) 



p>p l 
p2 



<5>(p) = N 



2mp — (p + 4m )arctanh(p/2m) 
2P 



(5) 



(6) 



(7) 



is the polarization tensor in the lowest, one loop, or- 
der. Here rj™ = diag(l, vf, vf), and p denotes the 
rescaled momenta pj = r^p^. N is the number of two- 
component fermion species, N = 4 for graphene. Note, 
that in chosen above representation of gamma matrices 
the parity-odd parts of the polarization tensor II are can- 
celed between contributions of various fermion species. 
If they are not, due to external magnetic field or for 
some other reason, this may be measured by studying the 
polarization rotation of light passing through suspended 
graphene 20 . It is also essential to notice that both po- 
larization operator II ([6]) and the diagram in ([3]) do not 
posses infra-red divergencies in the m = limit. 

To calculate the diagram ([3]) we only need to couple 
the kernel © to the photon propagator and integrate 
over the photon momenta. Symbolically, 



r = -Tr(riD), 



(8) 



where Tr is the functional trace, and D is the photon 
propagator. In the Feynman gauge the propagator 
is diagonal. Since IT- 7 by construction does not have com- 
ponents along the ir 3 direction, we are only interested in 
the Dij part of the propagator. According to ((4]), these 
components satisfy Dirichlet boundary conditions, and 
one can write 

Dij(x; y) = g. tj (D (x - y) - D (x - y R )) , (9) 

where Dq (x — y) is the standard (Feynman) propagator 
of a free massless scalar field in 3 + 1 dimensions, and 
the coordinate y R is reflected at the conductor surface, 
y 3 R = —y 3 . The full effective action T thus reads 

r= I J ^^y^ii x MDo(.x-y,0)-D (x-y,2a)], 

( 10 ) 

where both x and y lay on the surface of graphene. After 
making the Fourier transform in the directions parallel 
to the surfaces, one can write Dq as 



Do(x,y) 



d 3 p 



e iPi(^-y j ) D(p,x 3 -y 3 ). (11) 



For the Euclidean 3- momenta, i.e., after the Wick rota- 
tion p — > pe — (pi,p\ 1 p2), Pi — ipQ, its explicit form 
reads 



D( PE ,x 3 -y 3 



P\\ = \pe\- (12) 
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The term with Dq(x — y, 0) on the right hand side of (fT0|) 
does not depend on a and will be neglected (as it does not 
contribute to the Casimir force). The remaining terms 
in r are non-divergent. After the Fourier transformation 
and the Wick rotation we obtain 



£1 



r 

TS 



d 3 PE m( PE ) 



-2ap,, 



2tt) 3 

1 f d 3 p E «(pf - 
4 J (2^ 



P\\P\\ 



-2ap|| 



(13) 



where we expanded lK(ps) explicitly with help of 



III. LIFSHITZ FORMULA APPROACH 

One can also adopt another point of view on the 
Casimir effect for this system and consider an effective 
theory of the electromagnetic field described by the ac- 
tion Sm + <Scff subject to the conductor boundary condi- 
tions (|4]) at x 3 = 0. Away from the surfaces, the photons 
propagate freely. They are reflected at the surface of 
the conductor x 3 = 0. At the surface of graphene, the 
Maxwell equations receive a singular contribution 



dpF*" + 6{x 3 - a)W p A p = 



(14) 



following from S e g. Here we extended II to a 4 x 4 matrix 
by setting II 3 - 7 = IP 3 = 0. This singular contribution is 
equivalent to imposing the matching conditions 



.4 



A, 



fi\x 6 =a+0 — ^fi \x i =a — Oi 

(d 3 A fJi )\ x 3 =a+0 - (<9 3 A«)U 3 =a-o = 11^ A v \ x3=a (15) 

At this stage, one can forget the origin of IL^" and quan- 
tize the electromagnetic field subject to the conditions 



3 at x 3 



and to the conditions (0| at x 3 = 0. 
This can be done, at least at a somewhat formal level, 
and even the renormalization theory can be developed, 
e.g. along the lines of Refill. The Casimir energy den- 
sity can be defined as a sum over the eigenfrequencies 
which then are expressed trough the scattering data of 
the electromagnetic field. There are several versions of 
this procedure leading to different representations for the 
Casimir energy - for details see Ref A 

On the other hand, the Lifshitz approach 2 ^ relates the 
Casimir energy density for two parallel dielectric slabs to 
the corresponding dielectric permittivities taken at the 
imaginary frequency. In some later work o 23 i 24 the con- 
nection between two approaches was established and gen- 
eralization of the Lifshitz formula was presented. For the 
interaction between two plane parallel interfaces sepa- 
rated by the distance a and possessing arbitrary reflection 



'TE.TM' 



(2) 
r TE,TM 



of the TE and TM electro- 



magnetic modes on each of the surfaces it reads 



£l = 



d 3 p E 
16tt 3 



ln[(l- 



-2p|,a (1) (2)w, 



-2p|,ajl) (2) 



(16) 



The reflection coefficients are to be found from cor- 
responding boundary or matching conditions. For 
graphene with help of (fT5|) we obtain at the Euclidean 
momenta 



'te 



2pi| + a<I> 



'tm 



ap\ 



2pf, + apii$' 



while for the perfect conductor one has 



(2) 
r TE 



-1. 



(2) 



1. 



(17) 



(18) 



It is clear, that $ must be rotated to Euclidean momenta 
as well. We also note that the perfect conductor case is 
recovered from (|17|) in the formal limit $ — > oo. 

The Euclidean momenta representation of the Casimir 
energy like (fl"6 |) -([T8 |) . has several advantages. First of 
all, it takes into account contributions from possible sur- 
face plasmon modes not requiring to analyze such modes 
explicitly^. Secondly, it is straightforward to consider its 
limiting cases and also perform numerical evaluations as 
presented in the next section. 

One can show by a direct computation that the energy 
£ i , Eq. (|13p , coincides with the leading a 1 order in a per- 
turbative expansion of the Lifshitz formula (fl6 |l - (rr8|) . so 
that the two approaches are consistent. In fact, the Lif- 
shitz formula is the one-loop vacuum energy (one closed 
vacuum loop) in an effective theory corresponding to the 
action Sm + S e g. As we have explained above, the pres- 
ence of the singular part S e g is equivalent to the match- 
ing conditions (fT5| . Imposing these conditions, in turn, 
is equivalent to summing up the photon propagators with 
an arbitrary number of n m ™ insertions (or, with an arbi- 
trary number of the fermion loop insertions). Therefore, 
the Lifshitz approach corresponds to a partial summa- 
tion of diagrams of the QFT approach. We shall explain 
this correspondence in detail elsewhere. 



IV. RESULTS AND DISCUSSION 



The formulae (TIB")) and (fT6| - (|18p are suitable for the nu- 
merical and asymptotical evaluation. First we consider 
the large separation limit, a — > oo. Introducing dimen- 
sionless variables pe — * Pe& in IT6 l) -(fT8 l) we are able to 
expand the integrand in a power series in 1/a. Each term 
of such expansion is integrable and in the leading order 
we get for the energy 



aN 2 



J f 



96tt 2 



(19) 



We note that the energy is decreasing by one power of 
the separation a faster than for ideal conductors. Also we 
point out that this asymptotic expression is of the first 
order in the coupling constant a. 

In the limit of small separation, a — > 0, due to the 
structure of the function $(p^/a) the energy 



1 



a ^ 167T0 3 



h(a, N, vf) 



(20) 
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factorizes into a distance dependent part and a function 
h(a, N, Vp) independent of the separation a but contain- 
ing all powers of a; in the leading order it is 



Na 

h(a,N,v F ) = — — x 
lb 

ft 2 + V F 

x 1 H , arcsmh 



(21) 



0^ 




+ 0(a 2 ). 



Therefore we see that the asymptotic behavior of the 
Casimir energy in our model shows some surprising fea- 
tures being drastically different from that in the hydro- 
dynamic model 5 - 7 . At large separation it does not turn 
into the ideal conductor case, while at small ones this 
case is regained. Such behavior is counter-intuitive since 
the main contribution at small separations shall come 
from the high momenta for which one would expect the 
graphcnc film to become transparent. On the other hand, 
the behavior at large separation results from low momen- 
tum contributions for which the Dirac model is consid- 
ered to be well proved following directly from the elec- 
tronic structure of graphene. 

One can show that in the case of a massless (gapless) 
fermions the same distance dependence as for ideal con- 
ductors is retained for all distances due to lack of any 
dimensional parameters in this case. The magnitude of 
the energy is defined by the same asymptotic ([20]) as for 
the a — > case. 

Let us now turn to the numerical evaluation. It is 
convenient to normalize the results to the Casimir energy 
density 



Sr 



720 a 3 



(22) 



for two plane ideal conductors separated by the same 
distance a. The relative quantities £i/Sc ancl £l/£c are 
dimensionless and depend on a single dimensionless pa- 
rameter ma. To fix the scale, note that for m = O.leV 
(actual values of m are much smaller) ma — 1 corre- 
sponds to a = 1.97 micrometer. The results of calcula- 
tions are depicted at Fig. Q] For m — the normalized 
energies S^l/Sc are constant independent of a as ex- 
plained above. 

Thus, we can see that the magnitude of the considered 
Casimir interaction of graphene with a perfect conductor 
is rather small. Actual measurement of such weak forces 
is a challenging, but by no means hopeless, experimental 
problem 31 . Strong dependence on the mass parameter 
m at large separation is also a characteristic feature of 
the Casimir force. Getting an independent measurement 
of m may be very important for our understanding of 
the electronic properties of graphene. The mass of quasi- 
particles in graphene is, probably, very tiny. This im- 
proves the detectability of the Casimir interaction since 
the energy increases with decreasing m. 

As noted above, the Casimir energy £i calculated in 
the QFT approach coincides with the lowest order ex- 
pansion in the coupling constant a of the Casimir energy 




ma 



FIG. 1: The relative Casimir energy densities £i/£c (solid 
line) and £l/£c (dashed line) as functions of ma. A zoom of 
the small-distance region is placed in the upper-right corner. 



£l derived from the Lifshitz formula. This perturbative 
nature of the Casimir energy in graphene systems (which 
is just the other side of its smallness) makes the calcu- 
lations much easier and will probably simplify the anal- 
ysis of other geometries like, e.g., folded or corrugated 
graphene near a conducting surfaces. (An example of 
perturbative calculations with ^-potentials with the sup- 
port of a non-trivial geometry can be found in Ref^ 5 -). 
This is in contrast with some other physical effects in 
graphene which exhibit a strong coupling dependence. 

There are several factors that can be included in the 
considered model to make it more realistic. Among 
them the effects of non-zero temperature, corrugation 
of the free-standing graphene sample, presence of impu- 
rities and non- vanishing density of carriers. However, 
for clean enough graphene samples, similar to Refi2£, all 
mentioned factors can be considered as perturbations not 
changing the essentials of the Dirac model. In particu- 
lar, the corrugation was shown to maintain the massless 
fermionic nature of the quasi-particles^ 7 -, and thus can be 
treated perturbatively as in Refi^ 5 .. The effect of impu- 
rities can also be implemented into the theory as modi- 
fications to the Dirac operator which may be considered 
as perturbations at least in the sense of the operator the- 
ory. Indeed, while the impurities have a significant im- 
pact on the density of states^, and broaden the Landau 
levels in the Hall regime, their effects may be represented 
through adding a complex chemical potential, expressed 
through electronic self-energy^. Although in general the 
frequency dependence of these terms can significantly 
modify the properties of the model, in many cases the 
constant approximation of such broadening of Landau 
levels works well2^ and can be treated as corrections 3 ^. 

Thus, the method proposed above will therefore re- 
main valid upon inclusion of these effects, and probably 
even numerically the corrections will remain small. Ac- 
tual calculations will be the scope of our future work. 
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